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' We present an extensive comparative study of ground-state densities and pair distribution func- 

' tions for electrons confined in two-dimensional parabolic quantum dots over a broad range of cou- 

' pling strength and electron number. We first use spin-density-functional theory to determine spin 

, densities that are compared with Diffusion Monte Carlo (DMC) data. This accurate knowledge of 

f— I ' one-body properties is then used to construct and test a local approximation for the electron-pair 

' correlations. We find very satisfactory agreement between this local scheme and the available DMC 
data, and provide a detailed picture of two-body correlations in a coupling-strength regime preceding 

■ the formation of Wigner-like electron ordering. 
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I. INTRODUCTION 

^ . For a number of years there has been a growing interest in studying finite quantum systems under external confine- 
(— ! ■ ment, such as ultracold atomic or molecular gases inside magnetic or optical trapsi and electrons in metallic clusters^ 
or quantum dots^. The confinement introduces a new length scale and induces novel physical behaviors relative to 
the corresponding infinitely extended model system. In particular, in a quantum dot the properties of a homogeneous 
electron gas are profoundly modified by the emergence of effects that are commonly associated with electrons in 
atoms. A well-known example is the presence of a shell structure in the energy to add electrons to a quantum dof*. 
' With regard to the spatial structure of the electronic system, the analogue of two-dimensional (2D) Wigner crys- 
tallization has been shown in a path-integral Monte Carlo study^ to occur in two distinct stages inside a circularly 
O symmetric parabolic quantum dot. Radial ordering of the electrons into shells occurs first and is followed by ori- 
, O, . entational ordering through freezing of intershell rotations. Short-range order in the electronic structure at lower 
coupling strength is described by the pair distribution function g{r, r') giving the spin-averaged probability of finding 
two electrons at positions r and r'. Some properties of this function and of its extension to describe spin-resolved pair 
J> ' correlations have been reported for a circular quantum dot in a DMC study by Pederiva et al^. In the macroscopic 
limit g(r,r') reduces to a function of the relative distance |r — r'| of an electron pair and describes the liquid-like 
I short-range order in the homogenous electron gas. 
^ j ' The main purpose of this work is to present an approximate theoretical treatment of the electron-pair correlations 
( <~>) in quantum dots at such values of the coupling strength and of the electron number. We use spin-density-functional 
ly-j theory^ (SDFT) in local spin density approximations and test their accuracy against the available DMC data of 
, Ref. i6| for both one-body and two-body structural properties. In brief, the paper is organized as follows. Section Hll 
■ summarizes the SDFT procedure for the sake of completeness and in order to set the stage for the later sections of 
' the paper. Our results are presented and discussed in Section imi for the one-body spin densities and in Section HVI 
H for the pair correlation functions. An Appendix emphasizes the distinction between the short-range order of present 
I ' interest and the broken-symmetry states that are met in unrestricted Hartree-Fock (HF) calculations* on electrons in 
quantum dots. Finally, a brief summary of our main conclusions is given in Section El 
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II. THEORETICAL APPROACH 



H real-space Hamiltonian is 



We consider N interacting electrons of band mass m, confined in a strictly 2D parabolic quantum dot (QD). The 



+ ^E/^''-/ rfVvii(r)^l,(r')«(|r-r'|)Vi.'(r')V'.(r). (1) 

cr.cr' 

Here V'o-(r) and ip^iv) are Schrodinger field operators obeying canonical anticommutation relations, VJ;xt(r) = mujQr'^/2 
is the external confining potential and v{r) = e^/^nr) is the interparticle Coulomb potential, k being the dielectric 
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constant of the material. The Hamiltonian Q commutes with the z-component Sz of the total spin and therefore 5*2 
is a good quantum number. 

We introduce the spin density n„[r) — (-0j(r)i/)cr(r)), the density matrix p„„'{t,t') = {'ip'l{v)'ij}„i[r')), and the 
two-body density n^^j,(r,r') = (i/it (r)^^, (r')')/)^/ (r')V)^(r)) = n^(r)n„/ (r')^^^/ (r, r'), g„„,{r,r') being the spin-resolved 
pair distribution function (PDF). The calculations will be carried out at fixed N = + Ni and Sz = {N^ — Ni)/2, 
where — J (Prna-{r) is the total number of electrons with spin a i- 

Choosing as unit of length the harmonic-oscillator length £o = y/h/{mu!o) and as unit of energy the harmonic- 
oscillator quantum hioo, the QD Hamiltonian can be shown to be governed by the dimensionless parameter 

where = kY? j (me^) is the effective Bohr radius. The physical properties of the electron assembly are thus functions 
of the quantities A, iV, and Sz- We will choose parameters that are suitable for a 2D electron gas (EG) confined in 
a GaAs quantum well, i.e. m — 0.067 bare electron masses and k = 12.4. With this choice a*^ « 9.8 nm and the 
effective Hartree energy is e^/(Ka^) w 11.9 meV. 



A. Spin-density-functional theory in the Kohn-Sham scheme 

Within the Kohn-Sham (KS) version of SDFT the calculation of the equilibrium densities ^^.(r) is recast into the 
solution of a set of Schrodinger-like equations for the single-particle orbitals (r), 
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-f Kxt(r) + t;H(r; K]) + <,(r; {n„\) 



(3) 



Here i'H(r; \na\) — 'Ylia /d^r'?;(|r — v'\)n„(v'^ is the classical Hartree potential and v'^J^; [f^(T]) = 5E^t^(n„\| &n„(y) is 
the spin-dependent exchange-correlation (xc) potential defined as the functional derivative of the xc energy functional 
£'xc['^cr]- The approximation that we have employed for Wxc(r; \na\) is discussed below. 

The KS mapping guarantees that the equilibrium spin densities can be built from the KS orbitals in a single-particle 
fashion, 



.(r) = ICW 



(4) 



i (occ) 



where the sum runs over all occupied states. This equation also provides a self-consistent closure for the KS equations 
©. Once these equations have been solved, the ground-state energy E of the system is obtained from E — T^^- 
En [n^] + Vcxt [n^] + -E-xc [n^] , where Tg is the kinetic energy of the auxiliary noninteracting electron system. En [n^] is 
the Hartree potential energy, and Vcxti'T'o-] is the contribution from the external potential. 



B. Adiabatic connection and approximate xc potential 



An implicit expression for the xc energy functional, which highlights the importance of the PDF, is the adiabatic 
connection formula^. This reads 

EM = ^E/^'^/ rf'r't;(|r-r'|)n,(r)n,,(r')[W(r,r';K])-l], (5) 

where 

r'; KD = d^gi^lir, r'; [n,]) (6) 
Jo 

is the coupling-constant averaged PDF. Here g'^J^, is the PDF for a system with interactions v^{r) = e^f/(Kr) and 
fixed (^-independent) equilibrium densities ncr(r). This function depends on no-(r) as a result of the Hohenberg-Kohn 
theorem. 
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The local spin density approximationiS (LSDA) for E^dna] reads 

£^x"c'°^K] = / d'rniv) eL°'"(",C)L_(r).c-.c(r) ' (7) 

where e^°™{nX) is the xc energy per particle of a homogeneous 2D electron gas (EG) as a function of the total 
particle density n = X^o- '^<^ ^^'^ spin-polarization ( = X]<t '''^o-/'^- Accurate results for £^°'"(n, have been 

obtained in Quantum Monte Carlo simulations by Attaccalite et aL^- with special attention to its dependence on (. 

The adiabatic connection formula allows one to interpret the LSDA as an approximate choice for g^g-i (r,r'; [n^]). 
One readily obtains Eq. Q by approximating the xc energy density from the exact expression in Eq. ((SJ with 
ne^°'^{n, C) taken at the local density n{r) and at the local spin-polarization C(r)- The xc energy of the EG is given 
by 

s^rin, C) = ^ E / v{r) [g^fj}\r; n, () - 1] , (8) 

(7,(7' ' 

where g^°™(r; n, C) is the coupling-constant averaged PDF in the 2Z? EG. Thus within the LSDA the exact functional 
gcrcr' (r, r'; [rio-]) is approximated as 

g^^'(r,r'; [n^]) ~ .g;j°7'(|r - r'|; n, C)|„_^„(r)^(-_^(r) = W (r, i")Ilsda • (9) 

Analytical representations of accurate Quantum Monte Carlo data for the spin-averaged PDF of the homogeneous 
2_D EG are available in the literatures^ and provide a convenient input for our work. 

We recall at this point for later use that the exact PDF satisfies the so-called central sum rule 

y d\'n,{r') [g„,,n (r, r') - 1] = -6,,, 6,^,, , (10) 
in addition to the symmetry property g^ai (r, r') = g^ia{v' , r) and to the asymptotic result lim|r_r'|^-i-oo So-o-' (r, r') — 1. 

C. Fock-Darwin basis 

In the numerical solution of the KS equations we have adopted a standard procedure involving projection of Eq. 
on the Fock-Darwin (FD) basis corresponding to the (complete and orthonormal) set of eigenfunctions of the 2D 
isotropic harmonic oscillator. These are the product of the eigenstates of the angular momentum — —ihdg and 
of the radial functions Rn,Mif), ¥'n,Af(r) = (27r)~^/^ exp (iAf0)i?„^M('')- The quantum numbers n and M represent 
the number of nodes of Rn,M{f) and the angular momentum Mh carried by the state. The radial wave functions are 
expressed through the generalized Laguerre polynomials^. 

The projection of Eq. (0) onto the FD basis is straightforward. Decomposition of the KS orbitals, (/)f^(r) — 
Tlict C^'°'(^Q(r) where a stands for the pair {n^, Mq,}, leads to a matrix eigenvalue problem.^'* for the coefficients C^''^, 

E(«|HKs|/3)C^'^=e.,.Cr. (11) 

Here TYks is the effective KS Hamiltonian in Eq. Q. The equilibrium densities in Eq. |0J take the form 
"-ct(i") — X^a /3(Si[CQ'^]*C'^''^)'/'a('')'P/5(i"): which is uscd in evaluating the xc potential from Eq. Q. The Hartree 
term in Eq. © is expressed through (a|uH(r; [ricr])|/3) = Y.i,a Y.i,8[CY^*'~^l''^ Va^sp, where Va.ysp = / d'^vd'^r' v{\v - 
r'|)(/3* (r) (y9* (r') (/^^(r') (/7/3(r) are the two-body Coulomb matrix elements. Selection rules on the quantum numbers 
are hidden in V^^a/?: for instance, Va-ySfi vanishes unless the angular momentum is conserved in a scattering process, 
i.e. unless M^ + M-y = Ms + Mp. This is easily verified through an expansion of l/|r — r'| in cylindrical coordinatesi^. 

In practice, the sums over the FD basis elements must be truncated. The numerical calculations have used A^max = 
20 energy levels, which corresponds to (A^max + l)(A'max + 2)/2 = 231 single-particle states. Convergence of the 
self-consistent procedure has been achieved with a precision of at least 10"^ on the electron density. 
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III. TESTING THE LSDA FOR THE ONE-BODY DENSITIES 



A main aim of this section is to compare our LSDA results for the one-body density profiles with the state-of-the-art 
DMC data of Pederiva et al&. The comparison confirms the conclusions already drawn in Ref. and will give us 
confidence in the inputs to be used in our calculations of electron-pair correlations that will be reported in Sect. Hvl 
We also report numerical results obtained within the HF approximation^ and in some cases (iV = 6 and 30) we 
illustrate the role of electron-electron interactions by showing the density profiles for noninteracting electrons (A = 0). 

Our LSDA and HF calculations of na-{r) for circular 2D QD's refer to the cases N = 6,9,12,20 and 30. The 
confinement energy has been chosen as hwo = 3.32 meV, which corresponds to A = 1.89. A summary of our main 
results is shown in Figs. 11131 

It is immediately evident from Fig. ^ that the LSDA density profiles are in excellent agreement with the DMC 
data, except for = 12. We have no explanation for this specific discrepancy. In the case A^ = 9 the ground state 
at A = 1.89 is partially spin-polarized with Sz = 3/2, while in all other cases it is paramagnetic {Sz = 0). Figures[21 
and O show that in the spin-polarized case at A^ = 9 the agreement with the DMC data is excellent for both the total 
density profile n{r) and the local spin polarization C,{r). We have also checked that our LSDA results are not unduly 
sensitive to the input chosen for ex°™(n, Q). We have tested in this respect the earlier parametrization of ex°'°(n, C) 
given by Tanatar and Ceperle}»i& and found minor differences arising in the local spin polarization for A^ = 9, as is 
shown in the inset in Fig. |3| 

We also confirm that the HF is not a good approximation for the ground-state density profiles, especially for small 
values of A^ where the role of correlations is more important. The quality of the HF results appears to improve with 
increasing A^, as indicated by the case A^ = 30 in Fig.^ The local spin polarization for the case A^ = 9 in Fig.|21is 
also reasonably accounted for. A brief discussion of symmetry breaking in HF calculations of the one-body density is 
given in the Appendix. 



We have seen in Sect. Ill 51 how the PDF enters the adiabatic connection formula for a formal definition of the xc 
energy functional. The PDF directly describes the conditional probability density P(r', cr'|r, cr) — ncr'{r')gcra'{r,r') of 
finding an electron with spin a' at position r' when another electron with spin a is at position r. In the homogeneous 
2D electron fluid the increase of coupling strength with decreasing particle density towards a spin-polarized state 
and a triangular Wigner crystal is accompanied by strengthening short-range order in the electron-pair distribution. 
This is signaled by the emergence of a peak in g{r) at a relative distance r approaching the first-neighbor distance 
dwc = (%/3n/2)-i/2 in the crystalii*i&. 

Of course, the standard formulation of SDFT only gives access to the equilibrium one-body densities ^^(r). Several 

attempts have been madeii to build a generalized functional approach having the pair density n^^J, (r, r') as its basic 
variable, from which both n^ir) and gcrcr'(r, r') may be obtained. A practicable self-consistent procedure to calculate 
(7ctct'(''i been proposed by Davoudi et ali^, who extended to inhomogeneous fluids an approach originally used 

by Overhauser^iS to evaluate electron-pair correlations at contact. An Overhauser-type approach has also been set 
up2£ for calculating the angularly and center-of-mass averaged pair density, which suffices for evaluating the xc energy 
of an inhomogeneous electron system. 



In the present context, we examine an alternative approximate approach to the PDF, allowing relatively simple 
numerical calculations with results that will be compared with the DMC data of Pederiva et al& for a 2D QD. Our 
approach is inspired to the so-called average-density and weighted-density approximations, that have been proposed 
in the literature for the purpose of transcending the LSDA in the evaluation of the xc energy functional (sec Drcizlcr 
and Grossi and references therein). These approximations satisfy by construction the "central sum rule" in Eq. 110|l . 

In this so-called average-spin-density approximation (ASDA) the functional dependence of gaa'{i^, i"'; [n-a]) on na{r) 
is taken in the form 



IV. PAIR CORRELATIONS 



A. The average spin-density approximation for the PDF 



w(r,r';K])c. .g^r(|r-r'|;n,C) 



n— *n(r,r'),(^— >(^(r,r') 



gaa' {r,r')\ 



ASDA 



(12) 



where 



n(r,r') = [n(r) n(r')]/2 
C(r,r') = [C(r) + C(r')]/2 



(13) 
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(see also the work of Ebner et alm^ and of Yamashita and Ichiman*22i) . Contrary to the LSDA, the ASDA satisfies 
the symmetry property go-o-' (''j = da'ai^' t^)- It still implies some minor violations of the central sum rule: for 
instance, we have verified that it may lead to deviations from the requirement in Eq. 1)10(1 which become as large as 
5% for the spin-summed PDF in the bulk of a QD with = 9 and A = 1.89. 

In Fig. 2] we compare our ASDA and LSDA results for the spin-summed PDF, defined as 



^-^ n(r)n[r) 



a, a 



with the DMC data of Pederiva et al.^ on a QD with = 9 electrons and A — 1.89. The quantity being shown in 
Fig. 01is g{r,r' — 0), which depends only on the modulus r — |r| owing to the circular symmetry of the ground-state 
density. As already noted, the ASDA satisfies the symmetry property of g^-a'i^,^') and it does not matter whether 
one sets r or r' to zero (corresponding to the center of the QD). However this is not the case for the LSDA, and we 
have decided to show in Fig. ^ the choice that corresponds to setting r' = 0. 

It is seen in Fig. 01 that the ASDA and the LSDA give essentially the same results when the second electron is 
also close to the center of the QD, and are in good agreement with the DMC data. But the LSDA badly fails 
in describing long-range correlations, because it breaks down across the edge of the QD at r ~ 3£o where n{r) is 
rapidly dropping to zero. The ASDA is instead calculated at the average density n(r, 0), which tends smoothly to a 
constant across the QD edge. These behaviors can be emphasized by referring to local definitions of the density 
parameter as r^^°*(r) — [TTn{r)]^^^'^ /a*^ and r^^°*(r, r') = [7rn(r, r')]~^^^ /a^. As is shown in the inset in Fig.^ while 
r*^°*(r, 0) remains essentially constant on crossing the QD edge, r^^°*(r) increases in an exponential way heralding 
the breakdown of the LSDA. In practice, however, this breakdown is less serious than it may seem, since the two-body 
density is determined by the PDF multiplied by density factors. 

We proceed to present a broader view of the ASDA spin- averaged PDF for the same partially spin-polarized QD. 
Figure m shows the geometrical coordinates that will be used in the following figures. In Fig. we show a three- 
dimensional plot of g{r, r'), when both r = xx and r' — x' x lie on a line x passing through the center of the confining 
potential (see Fig. [S] left). The main features in Fig. are as follows: (i) the Pauli-Coulomb hole lying along the 
diagonal x = x'; (ii) the correlation-induced oscillations which are seen to lie along directions parallel to this diagonal, 
as are better seen in the contour plot; and (iii) the essentially flat asymptotic regions further out. We may remark that 
r^^°^{x,x') is a smooth and bounded function if at least one of the two coordinates lies in the bulk of the QD inside 
its edge. The calculation loses meaning when both coordinates are far outside the edge, so that n{x,x') is rapidly 
vanishing. In fact, the needed input on the PDF of the homogeneous 2D EG from Rcf. 12 is limited to electron 
densities corresponding to = (7rn)~^/^/a^ up to the value 40. The contour plot in Fig. |51shows as hatched areas 
these regions of inapplicability, located approximately at (|a;|, |a;'|) > 3^o- 

Before concluding this section we should comment on the spin-resolved pair correlations for the same QD. Unfor- 
tunately Ref. 12 does not provide analytical representations for the spin-resolved PDF of the homogeneous 2D EG 
at finite values of (. This has prevented us from building the correspondent spin-resolved ASDA PDF for partially 
spin-polarized QD's. However, in Fig. [3 we compare the parallel-spin PDF's in the HF approximation, defined as 

g..(r,r)|Hp-l „^(r)„,(r') ' ^^^> 

with the DMC data of Pederiva et al.^. The quantity being shown in this figure is ga<j{r,0)\j^p. We conclude that 
at the value of the coupling strength in Fig. [7| the parallel-spin HF PDF is already in fairly good agreement with the 
DMC results. On the other hand, the HF approximation completely misses antiparallel-spin electron-pair correlations 
by giving 5|x(r,r')lHF ^ 1- 



B. Evolution of short-range order towards Wigner-like order in a 2D QD 

The attainment of Wigner-like order in confined electronic system has been studied by a number of authors. In 
particular, Egger et al2^ have reported a crossover from Fermi liquid to "Wigner molecule" structure in a finite- 
temperature study of QD's containing up to A^ = 8 electrons by path-integral Monte Carlo simulation (PIMC). 
A later PIMC study by Filinov et al.^ has regarded electron clusters in QD's with different particle numbers at 
various temperatures and coupling strengths. For even values of N these authors took the electronic system in a 
paramagnetic state and predicted a "phase diagram" at zero temperature, which shows a transition to a radially 
ordered state followed at much higher coupling strength by a transition to an angularly ordered state. For instance, 
in the case A^ = 10 (that is the lowest value of the particle number in their study) the first transition occurs at A ~ 20 
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and the second at A ~ 2770. The good predictive value of the ASDA for electron-pair correlations allows us to inspect 
how the short-range order in a QD at weak coupling evolves with increasing A towards radial Wigncr-like ordering. 
We do this below for the case = 6 and the results are presented in Figs. I8I11I 

In calculating the one-body radial probability density 27rrn(r) we enforce circular symmetry and consider only 
spin states that are paramagnetic [Sz ~ 0) or ferromagnetic (Sz = 3). The ground state changes spin polarization 
with increasing coupling strength: for A = 1.89 and 3.54 the paramagnetic state is lower in energy, but lies above 
the ferromagnetic state for A = 6.35, 10, and 12. Figure |H1 shows the probability density for both states of spin 
polarization at the above values of the coupling strength. A shoulder and ultimately a marked minimum appear in 
the ferromagnetic state with increasing A: similar results have already been reported by Egger et alm^ and by Reimann 
et alm^. The electronic system acquires the so-called (1,5) structure consisting of one electron at the center of the 
trap and a surrounding ring of five electrons. We have checked that Jq""" 2TTrn{r)dr = 1 for A = 10 and 12, where 
''mill is the position of the minimum in the probability density, and found that the height A(A) of the probability 
density at Tmin vanishes for A ~ 14 (see the inset in Fig. 

We turn to present the ASDA results for the evolution of the radial and angular dependence of the spin-summed 
PDF with increasing A. Figure reports the function g{r, 0) for the ground state of the QD and shows that, whereas 
the paramagnetic ground state at weak coupling does not possess any pronounced radial structure, the ferromagnetic 
ground state at A = 10 and 12 exhibits a main first-neighbor peak followed by secondary structures. All these 
structures are in phase with structures in the local coupling strength r^^°*(r, 0), as is shown in the inset in Fig. |^ 

The growth of radial ordering with increasing A in the ferromagnetic ground state is even more clearly displayed 
by plotting the total conditioned probability density 2'!rr P(r\r' = 0) = 2'!rrn(r)g{r,0), which carries information on 
both the one-body density distribution and the radial electron-pair correlations. This function is shown in Fig. [TUI 
for the QD under discussion. The value of the coupling strength A ~ 14 at which the height of the minimum in 
this function vanishes represents within the present theory our estimate for the location of the transition to radial 
Wigner-like ordering in a parabolic QD containing six spin-polarized electrons. 

Finally, the angular dependence of the electron-pair correlations in the ground state is illustrated in Fig. at 
the radial distance rmax corresponding to the location of the absolute maximum in the probability density. The 
function that is being plotted at various values of A is g{9) — 5^°'"(2rmax sin {0/2); n, C) evaluated at n = »T.(rinax) and 
C — C('"max), with 9 being the angle between r and r' as shown in the right panel of Fig. [5| Of course, the mapping 
between 17(6') and the PDF of the homogeneous 2D EG is a consequence of the ASDA. Even at A = 12 the angular 
ordering of the electronic system in the QD is seen from Fig. 1111 to be still very much liquid-like. Starting in the 
ferromagnetic state from an electron at 6* = on a circle at r = rniax, we find an enhanced probability of having a 
first neighbor on each side of it and two additional structures further out on the circle, but there is no evidence for 
an ordered fivefold ring of electrons. Indeed, the positions of the peak structures in g(9) are far from corresponding 
to regular pentagonal angles as would be appropriate for an angularly ordered (1, 5) structure. 

As a final remark we notice that, while the angular distance 9p{X) from the first-neighbor peaks in Fig. ^2 decreases 
with increasing A, the preferred first-neighbor distance c?(A) = 2rinax(A) sin (0p(A)/2) along the circle at Tmax is 
increasing with A. This is shown in the inset in Fig. 1111 and is due to the increase in Tmax with increasing Coulomb 
repulsions. 

V. SUMMARY AND CONCLUSIONS 

In summary, the main original parts of this work have concerned the theory of the short-range order that may be 
met in electron assemblies confined inside 2D parabolic quantum dots in a weak-to-intermediate range of coupling 
strength. We have proposed a very practical scheme for the calculation of the pair distribution functions in these 
inhomogeneous electron systems and examined in great detail its predictions in two specific cases. For a partially 
spin-polarized system of nine electrons we have seen that the theory is able to quantitatively account for the available 
Diffusion Monte Carlo data on spin-averaged two-body correlations. We have added to this a panoramic view of the 
charge-charge correlations, that waits to be tested in further Monte Carlo studies. 

The second problem that we have examined in detail has been the state of spatial short-range order in the paramag- 
netic and ferromagnetic states of a system of six electrons as a function of the coupling strength parameter. Naturally 
enough, by being based on a density functional approach that takes input from the homogeneous 2D electron gas, 
our predictions parallel to some extent the well-known phase behavior of this macroscopic system. On increasing the 
coupling strength the ground-state of the quantum dot first changes from paramagnetic to ferromagnetic and then 
acquires radial order in coexistence with orientational liquid-like short-range order, over the range of coupling strength 
that we have considered. It would be important, we feel, to re-examine these correlation properties in the quantum 
dot with six electrons by exact-diagonalization methods. One could test in this way to what extent the predictions 
that take their start from the macroscopic electron gas are in accord with those that are based on a few-electrons 
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atomic viewpoint. 
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APPENDIX: SOME COMMENTS ON BROKEN SYMMETRY STATES 



An approximate treatment of a strongly correlated many-body problem can in principle lead to states with sponta- 
neously broken rotational symmetry (see for instance the review of Reimann and Mannineni and the discussion given 
by Ring and Schucl*2^). A well studied example in the area of QD's is the self-consistent spin-and-space unrestricted 
HP treatment of the one-body density, which has been shownA to break the rotational symmetry at relatively low 
values of the coupling strength. 

As an illustrative example we show in Fig. ^] the HF one-body density for a QD with = 6 electrons at A = 3.18 
in both the paramagnetic and the ferromagnetic case, in full agreement with the findings of Ref. 0. The state of order 
that these pictures suggest for the six-electron system is very different from the results that we have found from our 
calculations of the electron-pair correlations, as reported in Figs. 181111 
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FIG. 1: Density profile n{r) (in units of £q^) as a function of r/£o for a paramagnetic QD with = 6, 12, 20 and 30 electrons 
at A = 1.89. The results of the LSDA and of the HF are compared with the DMC data of Ref. 0. The dash-dotted lines are for 
noninteracting electrons. 
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FIG. 2; Density profile n(r) (in units of £q ^) as a function of r/£o for a partially spin-polarized QD with A*' = 9 electrons at 
A — 1.89. The symbols are as in Fig.Q 



11 




FIG. 3: Local spin polarization ("(r) as a function of r /la for a partially spin-polarized QD with N = 9 electrons at A = 1.89. 
The symbols are as in Fig. Q The inset shows the difference between ({r) calculated with two different parametrizations for 
e5j?'"(n,C) (see text). 
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12 3 4 

r/4 

FIG. 4: Spin-summed PDF g{r, 0) as a function of r /lo for a partially spin-polarized QD with A'^ = 9 at A = 1.89. The results of 
the LSD A and ASDA are compared with the DMC data of Ref.|a The inset shows r^^''*(r) (dash-dotted line) and r*^°*(r, 0) 
(solid line) as functions of r/^o (see text). 




FIG. 5: Left: definition of x and x' appearing on the horizontal axis of Fig. |S|(r = xx and r' = I'x). Right: definition of the 
azimuthal angle 9 giving the abscissa in Fig. 1111 
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FIG. 6: Spin-summed PDF g{x,x') as a function of x/lo and x' /£q for a partially spin-polarizcd QD with A'^ = 9 electrons at 
A = 1.89. The bottom panel shows a contour plot of g(x,x'): the thin solid lines limit the bulk central region (|a;|, \x'\) < '3£o 
which is shown in the top panel. 




FIG. 7: Parallel-spin PDF gcra{r,0) as a function of r/£o for a partially spin-polarized QD with = 9 electrons at A = 1.89. 
The results of the Hartree-Fock approximation are compared with the DMC data of Ref. 0. 
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FIG. 8: Probability density (in units of ig^) as a function of r/£o for a QD with N = 6 electrons at varying A: profiles for the 
paramagnetic state (top) and for the ferromagnetic state (bottom). The inset shows the height A of the minimum (in units of 
£q^) as a function of A. 
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FIG. 10; Total conditioned probability density 2nr P(r\r' = 0) as a function of r/£o for the ground state of a QD with A'^ = 6 
electrons. The curves for A — 1.89 and 3.54 refer to the paramagnetic state, while those for A = 10 and 12 refer to the 
ferromagnetic state. 
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FIG. 11: Azimuthal plot of the spin-summed PDF as a function of the angle 6 defined in Fig. |51 for the ground state of a QD 
with A'' = 6 electrons. The reference electron is located at ^ = 0. The curves for A = 1.89 and 3.54 refer to the paramagnetic 
state, while those for A = 10 and 12 refer to the ferromagnetic state. The inset shows the preferred first-neighbor distance d(A) 
along the circle at r = rmax, in units of a*g. 
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FIG. 12: Rotationally broken HF one-body density (in units of £q'^) for a QD with N = 6 electrons at A = 3.18, in the 
paramagnetic (upper panel) and ferromagnetic (lower panel) case. At this value of A the HF approximation predicts the 
ferromagnetic state to lie at lower energy. 



